Triplet superconductivity in quasi one-dimensional systems. 
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We study a Hubbard hamiltonian, including a quite general nearest-neighbor interaction, 
parametrized by repulsion V, exchange interactions J Z ,J±, bond-charge interaction X and hop- 
ping of pairs W. The case of correlated hopping, in which the hopping between nearest neighbors 
depends upon the occupation of the two sites involved, is also described by the model for sufficiently 
weak interactions. We study the model in one dimension with usual continuum-limit field theory 
techniques, and determine the phase diagram. For arbitrary filling, we find a very simple necessary 
condition for the existence of dominant triplet superconducting correlations at large distance in the 
spin SU(2) symmetric case: 4V + J < 0. In the correlated hopping model, the three-body inter- 
action should be negative for positive V . We also compare the predictions of this weak-coupling 
treatment with numerical exact results for the correlated-hopping model obtained by diagonalizing 
small chains, and using novel techniques to determine the opening of the spin gap. 



I. INTRODUCTION 

In the last years, the study of extensions of the usual 
Hubbard model has, among others, two important mo- 
tivations. First, to explain many feaures of the rich 
diagram observed in the quasi-one dimensional organic 
compounds (TMTSF) 2 X (where TMTSF means tetram- 
ethiltetraselenafulvalcnc and X represents, PFg, CIO4 
or other complex), other interactions must be consid- 
ered in the model hanulipnian, in addition to the on- 
site Coulomb repulsionlla. Second, consensus increases 
about the fact that the usual Hubbard hamiltoniarfi, does 
not seem to define the minimal model for the super- 
conducting cuprates, while additional terms, might ex- 
plain thcj-physics of the superconducting phase of these 
materials! The proximity between the spin density wave 
(SDW) and the superconducting phases observed in the 
phase diagrams of both kind of systems .has been many 
times pointed out as a remarkable facfUm. The symme- 
try of the order parameter is, however, different for both 
systems. While it has n been observed to be d x 2_ y 2 for 
the case of the cupratesQ, experimental evidence suggests 
that the observed superconductivity in (TMTSF) 2 C104, 
and (TMTSF^PFg underpressures, as well as in the lay- 
ered compound Sr^RuOja is of triplet p-wave character. 
In this context, the extended Hubbard model with corre- 
lated hopping, is particularly interesting, as it seems to 
provide a good scenario for the occurence of both kind 
of superconducting instabilities, as well as the SDW and 
charge density wave (CDW) ones, dependma-oiuthe val- 
ues of the parameters and the filling factoii3li!J~lL3. 

The model is defined by the hamiltonian 

H= ( c h c jo + h - c -) {tAA{l ~ n is )(l - n js ) 

+ tsB n is n ]S + t AB [n lS (l - n ]S ) + n js {l - n is )]} 
+ U y^TtitTtjj, + V ^ n i n o ■ C 1 ) 

i <ij> 



where < ij > denotes nearest-neighbor sites. This 
model contains the one which— Hirsch proposed to give 
rise to hole superconductivity!^. In the context of the 
superconducting cuprates, Eq. (jl|) has been obtained as 
a one-band effective hamiltonian when a low energy re- 
duction of the extended three-band Hubbard hamiltonian 
is performed. For realistic parameters of the three-band 
model, calculated hopping parameters-satisfy the follow- 
ing relation: tAB > {tAA + tsB)/2EfrtZl. In particular, in 
the limit in which the Cu-0 hopping t p d is much smaller 
than the charge-transfer energies, tABr^ of the order of 
tpd while tAA and tsB are of order tp d £3. In two dimen- 
sions (2D), mean-field calculations, including the effect 
of spin fluctuations, support the existence of a supercon- 
ducting phase with d x i_ y i symmetry at finite doping, 
which is stabilized with the additioruof a next-nearest 
neighbor hopping in the hamiltonianla. While a realis- 
tic effective Hamiltonian for (TMTSF^X has not been 
constructed so far, Eq. (1) contains the main terms, if 
only the ground state with zero, one, and two particles 
in the singlet sector, of an adequately chosen cell are 
retainedtjH 

Collecting the correlated hopping terms can in one- 
two- and three-body terms, the hamiltonian reads 



i <ij> 

+ ( c lff c j'ff + H.c.)[-t + ta(«i<T + n ja) 



<ij>o 

+ t 3 n la n ja }}, 



(2) 



where t = t A A, h = t A A~t A B and t 3 = Hab ^aa ^bb- 
While two-body interactions are usual in many-particle 
problems, three-body interactions are more rare and in- 
troduce additional complications in most of the usual an- 
alytical many-body treatments. The appropriate mean- 
field reduction of the three-body term with parameter 
t 3 to effective two-body ones can be performed using 
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Wick's theorem and neglecting the resulting normal or- 
dered three body term, withpfckp-s/acuum representing the 
optimum Slater determinantll3'E2l. In ID, we have verified 
that in the continuum limit (representing the fermionic 
fields in terms of bosonic ones) , this procedure is equiv- 
alent to perform operator product expansions in the re- 
sulting haajiiltonian, keeping only relevant and marginal 
operatorsEj. In other words, both approaches are equiv- 
alent in the weak-coupling limit. The ensuing two-body 
Hamiltonian readsEZI: 

H eff = U^riitmi +V eff Y n i n j 

i <ij> 

+ (4 c ^ + H - c -)[* e// + A(n lff + n ja )\ 

<ij><7 



with 



= t-t 3 (3 T 2 -p 2 ), V eff = V + t 3 r, 

W = 2t 3 r, A = t 2 +pt s , J = 4t 3 r, (4) 

where r = (cl a Cj a ) and p — (ni a ) = n/2, with n being 
the number of particles per site. In ID and the weak 
coupling limit : 
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(5) 



The hamiltonian (g), with arbitrary interactions 
t e ff,A,W,J,U,V e ff, defines the most general model 
with nearest neighbor two-body interations, which con- 
serves charge and spin SU(2) symmetry. If in addition, 
the parameters are related by Eq. (Q) with V = 0, 
p = 1/2 and arbitrary r, the model also has pseudospin 
SU(2) symmetryEil. In the general case, we will also in- 
clude the possibility of anisotropic exchange (J z in one 
direction, J± in the other two), breaking spin SU(2) 
symmetry. Eq. (||) contains all the contributions up 
to nearest-neighbors of the Coul omb , interaction when 
written in the tight binding basisE3H. In this case, for 
weak screening of the interatomic repulsion, the rela- 
tion U > V eff >A>W~J>0 has been derived, 
while smaUe j . va lues of V are expected if the screening 
is efficicntHSEi Arbitrary values of the different inter- 
actions could, however, be expected when dealing with 
effective models, derived from some multiband model, as 
it is the case of Eq. (|j) . 

Several specific cases of the model Eq. (|^) have 
been studi£d_]^fbre using continuum limit field theory 
(CLFT^MEJ 2 ! In particular, the chain described by 
the correlated hoppiog model Eq. ([!]) was analized re- 
cently at half fillingO. However, no definite conclusions 



regarding possible dominance of superconducting corre- 
lations at large distances were obtained, and the exten- 
sion to other fillings remains open. An accurate phase 
diagram has also been obtained numerically using topo- 
logical transitionstl Other works an.models similar to 
Eqs. (0) and (§) are cited in Refs.BB. 

In this paper, we study the phase diagram of H e j f (j^) 

with generic parameters (t e ^ , t% , W, J 2 , J±, V e ^), in 
the weak-coupling regime, using the CLFT. The bound- 
aries of the region with dominant triplet superconducting 
correlations at large distances are given by simple ana- 
lytical expressions. This region will be denoted TS phase 
in the following. As a further step, we obtain the phase 
diagram of the correlated hamiltonian Eq. (Q), for weak 
and strong interactions, by numerical diagonalization of 
finite rings. The opening of the spin gap is detected ac- 
curately using a novel method based on results of confor- 
mal field theory and renormalization groupES, which in 
turn is equivalent to a topological transiton which cor- 
responds to a jump in a Berry phaseEj. The results of 
both approaches are compared and the conditions for the 
existence of the TS phase are discussed. We also analyze 
the other phases of the model, and discuss within which 
region of parameters, a superconducting instability with 
d x 2_ y 2 symmetry could be expected in 2D. Section II 
describes the weak coupling results. Results of the ex- 
act diagonalization of Eq. (|l]) and comparison with the 
CLFT are presented in section III. Section IV contains a 
discussion. 



II. CONTINUUM LIMIT FIELD THEORY. 



The CLFT, also called g-ology, is a weak coupling ap- 
proach. The whole procedu re . has be en explained in de- 
tail in many contributionsB'olialJJSEil. We, thus, present 
here only a brief explanation of the steps followed. The 
basic assumption is that the interactions are small, in 
comparisson with the Fermi energy. The non-interating 
energy dispersion relation is linearized around the two 
Fermi points and the interactions in the momentum space 
are expressed in terms of four different scattering pro- 
cesses, which are labeled by coupling constants g^ii (gn_) 
if they involve the same (opposite) spins. 

To describe the low-energy physics of the problem, the 
fermion operators are decomposed as 

c ja -> exp (ik F Rj)^ + ^(j) + exp (-ik F Rj)^ -, a (j), (6) 

where ^ r(J (j), r describe left- and right-moving 
fermions. In the continuum limit a — + 0, L — > oo, 
with aL finite, being a the lattice constant and L the 
number of sites of the lattice, these operators scale as 
^r.a(j) — * yfa^r.aix = ja). The hamiltonian (Q) can be 
written in terms of the continuunxijelds ty rt(7 (x). 

Using the same notation as VoitOu and neglecting irrel- 
evant operators, we obtain the following coefficients for 
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the different scattering processes of the hamiltonian Eq. 



Sill 

92\\ 

9i± 

92X = 

53_L 
34_L 



(2V eff + ^f)cos(2fc F a), 



J, 



54|| = V eSS + =f, 33|| = 

U + (2V ef / - y ) cos(2fc f a) - J± 

-2VF + 8Acos(fc F a) 

E. + yeff -h- 
2 4 

W + 4A cos(fc F a) 



Jz 



— - — cos(2fc F a) 



6 nA (U - 2V e " + ^ + J 



2 4 
4A cos(fcpa), 



2W0, 
VFcos(2fc F a) 



(7) 



where <$ ni i in the g%± (Umklapp) coupling contant indi- 
cate that this scattering process is active only at half- 
filling. The fermionic fields \P ri<T (x) can be written 
in terms of bosonic fields <j) p , </v, where p, (a) de- 
notes charge (spin) degrees of freedom, by recourse to 
a bosonization identityB and the hamiltonian can be ex- 
presed as: 



H eff = H„ 



Ha, 



(8) 



where H v ,(v = p, a) is a sine-Gordon hamiltonian: 

H v = v„ j dx{hjll{x) + {d x ^f\ 

+ cos(V 8nK u (p v )}, 
or 

where n„ is the moment conjugate to <p u , and 



(9) 



I', 



Vf 



< = «/K) 2 -(J) ! - '», 

.94|| + 34± 



TV 



33± 
V F = W F 



,7m, 



?; F = 2t e //sin(— ), K v = 



9i± 

27T ' 

34 1 1 - ff4i_ 
7T 

/ 2irv F + g v 



2ttv f 



9» 



(10) 



with 

9 p = 23i|| 



92\\-92±, g a 



251 1 1 - 92H +92J-- (11) 



The physics of the sine- Gordon _hapiltonian is well 
known from renormalization groupBEffil At half-filling, 
for |<73i_| < 9p the charge sector renormalizes to the 
Tomonaga-Luttinger fixed point, where charge excita- 
tions are gapless. Away form half-filling, charge ex- 
citations are gapless, except, for commensurate fillings 
n = p/q and strong repulsive interactions^ the system 
is insulating if K p < 1/q 2 . For example at quarter filling 



(n = 1/2), the critical value of K p is In terms 

of the parameters of Eq. (||), the condition for Luttingcr 
liquid behavior at half filling reduces to: 

2U + W ef f + 2 J x + J z < and 2V ef f < W) (12) 

For the particular case of the extended Hubbard model 
with correlated hopping (||), using Eqs. (Q),@ and (Q), 
this condition reads 

(U + 2V + 8t 3 /n < and V < 0) (13) 

in agreement with Ref.El. In particular, for V > the 
system has always a charge gap at half filling. 

Out of half filling, the hamiltonian H p (see Eqs. (^), 
(|To|), and (^)), reduces to a gaussian model (except 
for higher order Umklapp processes which are relevant 



apcnsurate fillings and K p < 1/4, as mentioned 



363). In this case, K p is not renormalized and if 



for co 
above 

K p > 1 superconducting correlations dominate at large 
distances. From Eq. (|l0[) , this happens when g p > 0. 
At half filling, K p is renormalized, but the initial value 
(Eq. (JTo|) ) should be larger than one in order to reach a 
final value larger than one at the gaussian point in the 
renormalization-group procedure. Thus, from Eqs. (|^), 
(|i"o|), (0) and (|il|), a necessary condition for the domi- 
nance of superconducting correlations at large distances 
is obtained: 



-U - W eff + (2V eff + y + J x ) cos(Trn) 
-2W-8Acos(7rn/2) > 0. 



(14) 



This condition is also sufficient away from half filling 
(n 1). The character of the dominant superconducting 
correlations at large distances (singlet or triplet) and si- 
multaneosly the opening of a spin gap is determined by 
H a (Eq. ^) and the flow of its paramenters under renor- 
malization. In the spin SU(2) invariant case (J z = J±), 
fja = 9i±- For negative g a a spin gap opens, the triplet 
superconducting (TS) correlations functions (CF) decay 
exponentially, while the singlet superconducting (SS) CF 
decay as d~ 1 / Kp with distance d. Instead, for positive 
g u , K a renormalizes to 1, there is no spin gap, the TS 
CF decay as d~ 1 ~ 1 / Kp In 1 / 2 d and dominate over the SS 
CF, which decay as d-^ 1 ^'" ln~ 3/2 J. Thus, in the 
Tomonaga-Luttinger liquid phase with J z = J± and 



K p > 1, using Eqs. (E 



1|), one sees that TS CF 



dominate at large distances if and ony if: 



J, 



U+ (2V e}} - |)cos(tto) - J 
-2VF + 8Acos(7m/2) > 0. 



(15) 



If the model is not spin SU(2) invariant, TS CF dominate, 
decaying as d~ x l Kf, ~ x l Ka when K p > 1 and g a > \gi±\, 
which in turn implies that the renormalized K a > 1. 

Adding Eqs. (14) and ( p^| ) a very simple necessary 
condition for the existence of the TS phase in the model 
Eq. (0) with J z = Ji is obtained: 
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W eff + J < 0. (16) 

This result is consistent with theoretical analysis which 
relate triplet superconductivity with ferromagnetismEil. 

For the model with correlated hoppin g in the form of 
Eq. (H), using Eqs. (||||) the conditions ( |l^ , ^5|) take the 
form: 

g 

— U + 2y(cos(7m) — 2) H 1% sin(7m/2) cos(-7rn) 

-8(h + tan/2) cos(7rn/2) > 0, (17) 

in order that superconducting CF dominate at large dis- 
tances and: 

U + 21/cos(7rn) i 3 sin(7rn/2) 

+8(t 2 + t 3 n/2) cos(7rn/2) > 0. (18) 

for the region in which the spin gap is closed. Adding 
both conditions leads to: 

V + -t 3 sin(7m/2) < 0, (19) 

as a necessary condition for the model to have a TS phase. 
In the extended Hubbard model (t 2 — tz = 0), Eqs. fll7| ) 
and Eqs. (|lj) imply that a TS phase can only exist only 
for 2/3 < n < 4/3. 

In addition to the SS and TS CF, the phase diagram 
at half filling is determined by the CF at large distances 
of the following order parameters for charge density wave 
(CDW), spin density wave (SDW), bond ordering W&tf 
(BOW) and spin bond ordering wave (SBOW) orderBO: 

Ocdw = y^(-l)*nj g 



~ cos{^2-KKp^p)cos{-s/2-KK a (j) a ), 

o S Dw = ^2(-^yvn ia 



~ sm( v / 2irKp4> p )sm(y / 2TrK (7 (j) rT ), 
Obow = ^^{c\+is c is + H.c.) 

i<7 

~ sin( *j2inr p <j)p ) cos ( \p2mK~ a ^a) , 
Osbow = ^2 a(c\ +ls c lS + H.c.) 

~ cos(^2TrK p (f> p ) sin(V27rif CT ^), (20) 

The opening of the spin gap is accompanied by the 
ordering of cj) s with expectation value (4> s ) = 0. When the 
charge gap is closed, except for the particular case g p — 
|i?3_l| for which K p renormalices to 1, g p renormalices to 
a positive value and K p > 1 at the fixed point. Then, as 
discussed above, TS (SS) CF dominate at large distances 
if the spin gap is closed (open). Instead, when the charge 
gap opens (|<?3x| > 9p) > 4>p orders with expectation value 
(<j> p ) = {{(j)p) = yJit/i&Kp) ) for negative (positive) 



g-3±. If the spin gap is positive, this implies CDW (BOW) 
order, while if the spin gap iSrGkssed the dominant CF are 
are the SBOW (SDW) onesB'cHl3. These considerations 
lead to the phase diagram shown in Fig. 1 for the general 
model Eq. ([|) at half filling in the isotropic case J± = 
J z = J. For W — J — 0, the model rcdufie&-to the 
extended Hubbard model studied previouslyQcJH, and 
no BOW or SBOW phases appear. For large positive 
U the leading power-law decay of the SDW, CDW and 
BOW CF is 1/d, but the logarithmic correction of the 
former (In 1 / 2 (24, makes the SDW the dominant CF at 
large distancea-l. The combination of parameters which 
control the existence of the BOW and SBOW phases is 
2W + J. The BOW (BSDW) CF dominate in a certain 
region of U and V if and only if 2W+J > (2W+J < 0), 
particularly for large U near the line V = U/2 — W— J/4 
at which the spin gap opens. 
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FIG. 1. Phase diagram of the general model Eq.(^) at 
half filling (n = 1) for J ± = J z = J. (a) 2W + J < 0, (b) 
2W + J > 0. 



III. NUMERICAL RESULTS 

In this section we calculate the phase diagram of the 
correlated hopping model Eq. ([!]) by exact diagonaliza- 
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tion of finite rings, and compare the results with those 
of the previous section. We concentrate on the electros-, 
hole symmetric case and set tAA = ^bb = L Numericalt2l 
and CLFTLi results for one particle per site (n = 1) have 
been done recently. We have concentrated on other two 
densities: n = 1/2 (quarter filling) and n = 2/3. The 
boundary of the superconducting phase was determined 
from the equation K p = 1 
expression: 



enough. However, these points are also consistent with 
the CLFT results. 



with K p calculated from the 



yJ™D p /2, 



(21) 



where the Drude weigth D p and the compressibility k 
were obtained extrapolating to the thermodynamic limit 
numerical results (obtained in the usual waye3) using a 
polynomial in l/L, where L is the length of the ring. For 
n = 1/2, we have used rings with L = 8, 12, and 16, and 
for n = 2/3 the lengths used were L = 6 and 12. We have 
also calculated K p using other two expressions which ia.- 
volve v p , the central charge and charge and spin gapsEj, 
to check for consistency and finite-size effects. The latter 
are in general very small, except for specific cases men- 
tioned below. The opening of the spin gap was deter- 
mined from the crossing of the levels of lowest energy in 
the sectors with total spin S = and S = 1 for periodic 
(antiperiodic) boundary conditions if N/2 is evenJodd), 
where N is the number of particles in the systemEa. This 
method is based on results of conformal field theory and 
renormalization group, which show in addition that at 
the crossing point, the finite size corrections of these ex- 
citation energies per site go as 1/L 2 , without logarithmic 
corrections. This allows a very accurate determination 
of the parameters for which the spin gap opens using 
finite-size scaling. This level crossing has also a topologi- 
cal significance, since at this point, the spin Berry phase, 
which can only take two, values: or 7r (mod 2tt) jumps 
in (it) at the transitioned. 

In Fig. 2, we show the evolution with doping of the 
phase diagram for V = 0. As soon as the density de- 
creases below n = 1, the Umklapp processes become ir- 
relevant and according to the CLFT, the degeneracy of 
TS and SBQi5pCF on one phase and the SS and CDW 
on the otherE2rll3 is broken in favor of the supercondunt- 
ing CF, since K p becomes larger than one. As doping 
increases, the SDW phase advances over the TS phase 
and at quarter filling (n = 1/2) the TS phase lies en- 
tirely within the region of negative U. Instead the SS 
phase advances rapidly over the BOW (or dimer ordered) 
phase at half filling. The agreement of the predictions of 
the singlet-triplet level crossing method for detecting the 
opening of the spin gap with the CLFT results is excel- 
lent at weak coupling (\tAB ~ 1| < 0-1), what confirms 
the accuracy of the method. The numerical results for 
the points at which K p = 1 also agree very well with 
the weak coupling results for n = 1/2. For n = 2/3, for 
which only two points (L = 6 and L = 12) were used in 
the finite-size scaling, and K p varies very slowly with the 
parameters near the non-interacting limit, the numeri- 
cal points in the weak coupling limit are not accurate 
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FIG. 2. Phase diagram of the correlated hopping model 
Eq.(|l|) for V = and several densities indicated at the top 
left of each figure. Solid squares indicate values of U and tAB 
for which K p = 1. The straight line K p — 1 according to the 
CLFT (Eq. (0) with t 2 = 1 - tAB, tz = -2t 2 ) is shown dot 
dashed. Solid circles are points at which the spin gap closes. 
The corresponding results according to the CLFT (Eq. (^)) 
are represented by the full straight line. 

It is remarkable that for n — 1, the CLFT results are 
quantitatively valid even at intermediate coupling. How- 
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ever, out of half filling and weak coupling, the CLFT 
overestimates the region in which the spin gap vanishes, 
and clearly underestimate the extension of the supercon- 
ducting phases. 
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FIG. 3. Same as Fig. 2 for U = and n = 2/3. The 
Dpen circles are points which correspond to the singlet-triplet 
crossing for L — 12 and deviate substantially from the ex- 
trapolated results that indicate the opening of a spin gap. 

In Fig. 3 we show the effect of V on the phase diagram 
for n — 2/3. We have chosen this density because a TS 
phase exists for positive U if tAB < 1- For tAB < 0.4 the 
results are affected by large finite size effects, perhaps be- 
cause of the proximity of phase separation (PS), and we 
were unable to obtain reliable results. For 0.4 < tAB < F 
increasing V, from zero, the line of K p = 1 is crossed 
first and then at larger V a spin gap opens, in quali- 
tative agreement with the CLFT results. However, we 
obtain that the spin gap opens also when V is decreased 
taking negative values. This crossing is out of the reach 
of the CLFT. As V is further decreased one expects PS. 
For tAB < 0.9 our results for the opening of the spin 
gap at negative V are affected by large finite size effects 
(reflected by the difference between the values of V in- 
dicated with solid and open circles in Fig. 3). This is 
probably caused by the proximity to. PS. The numerical 
investigation of PS is very delicateO and is beyond the 
scope of this work. 



IV. DISCUSSION 

We have studied the phase diagram of a Hubbard 
model Eq. (|j|) in the weak coupling limit, generalizing 
previous studies which use the continuum limit field the- 
ory. The model includs the most general form of nearest- 
neighbor two-body interactions which conserves charge 
and spin. The phase diagram is very rich, and six differ- 
ent phases can appear, according to the dominant corre- 



lation functions at large distances (see Fig. 1). In partic- 
ular, in the isotropic case, if there is a nearest-neighbor 
attraction (V < 0) or ferromagnetic exchange ( J < 0) a 
phase with dominant triplet superconducting correlations 
(TS phase) can exist. Specifically 4V + J < is a neces- 
sary condition for the existence of the TS phase. Work on 
weakly coupled chains using perturbation theory suggests 
that small interchain hopping stabilizes a 3D long-range 
order with finite critical temperature, which corresponds 
to the dominant correlations at large distances in the 
purely ID casec3. Thus, we expect that our results can 
be applied to real quasi one-dimensional materials, which 
can be described by an effective Hamiltonian like Eq. (jl]) 
or (||) . Since in real materials V is expected to be repul- 
sive, an efficient screening of the interatomic repulsions 
and effective ferromagnetic exchange (like that present in 
one dimensional cuprates containing edge-sharing CUO4 
unitsES) would be necessary conditions for the existence 
of triplet superconductivity. It has been proposed that 
in some ideal limit, ferromagnetism and triplet supercon- 
ductivity might be related by symmetry operations of the 
group SO(5)li 

For weak coupling, the general model Eq. (||), with 
parameters satisfying the relations Eq. ([!]) describes the 
correlated hopping model Eq. (|l|) . For this model to dis- 
play a TS phase in the weak coupling limit, it is neces- 
sary that V+(2/ir)(2t A B-tAA-tBB) sin(nn/2) < 0. We 
have studied this model beyond the weak coupling regime 
by numerical diagonalization of finite rings. In spite of 
the fact that the size of the studied systems is small, 
the results agree very well with those obtained with the 
field theory in the weak coupling regime. In particular, 
the accuracy of determining the gap by the— method of 
the crossing of singlet and triplet excitationsE3 (which in 
turn is equivalent to a topological transition in the spin 
Berry phasec3) is confirmed. For values of the correlated 
hopping which are outside the reach of the field theory, 
the regions of dominating superconducting correlations 
at large distances extend beyond the predictions of the 
weak coupling treatment. 

There are other physical phenomena which are outside 
the scope of the field theoretical treatment we followed. 
One of them is the opening of the spin gap for tAB < 1 
and negative V found in our numerical calculations (see 
Fig. 3). Another one is phase separation. In addition, for 
very small t&B (tAB = 0.2), there is numerical evidence of 
peaks at incommensurate wave vectors in char.ge-charge 
and spin-spin correlation functions at half fillingcil. These 
can be qualitatively understood using the formalism of 
the exact solution for t A B — Oej: roughly, the doubly oc- 
cupied sites are represented by effective bosons and the 
singly occupied sites by effective fermions. The Fermi 
wave vector of these effective fermions depends on U and 
hence, it is in general different from the non-interacting 
Fermi wave vector. Unfortunately, a quantitative analyt- 
ical calculation of these correlations functions for small 
tAB 7^ 0, is very difficult due to the huge degeneracy for 
tAB = 0. On the level of the field theory, one might 
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speculate that irrelevant operators, which we have ne- 
glected, become important at large couplings in an ad- 
equate renormalization group treatment and lead to in- 
commensurations. In addition to the TS phase, the re- 
gion of singlet superconductivity (SS) for tAB > 1 and 
U > (see Fig. 2) is particularly interesting. In the case 
<3 = 2tAB — tAA — tsB = 0, SS has been proposed and 
found in a mean-held treatment by Hirsch as a model for 
hole superconductivity^, and coqipmeiLby other numer- 
ical and analytical calculationsEj'OEOO. However, the 
case we studied numerically in section III is electron-hole 
symmetric and t% plays an important role. It is also essef 
tial in 2D to give rise to <i a .2_ a 2-wave superconductivity 
In the limit in which tAB is much larger than all other 
energy scales of the system, at half filling, a reasonable 
approximation to the ground state is obtained splitting 
the ring into consecutive dimers, and solving the Hamil- 
tonian in each dimer. Including the hopping between 
dimers in second-order perturbation theory leads an en- 
ergy which is above the energy calculated with density- 
matrix renormalization group by 1.6/u!j. This suggests 
a picture in which the system is composed of dimers, 
which behave as hard core bosons, being frozen at half 
filling (leading to a dimerized phase with long range or- 
der, BOW in Fig. (2)) but which aquire mobility out of 
half filling, giving rise to dominant SS correlations. If 
this image can be extended to 2D, we expect some kind 
of short range resonance-valence bond ground state at 
half filling, since dimers can be ordered in many different 
ways, and SS of s- or d x 2_ y 2-wave symmetry would natu- 
rally arise as the system is doped. The favored symmetry 
depends on the topology of the Fermi surface. Accord- 
ing to mean-field calculations and including spin fluctua- 
tions, the c? 2 ,2_j ( 2-wave symmetry is favored for moderate 
doping if a negativej-pext-nearest-neighbor hopping is in- 
cluded in the modeB. We are presently investigating this 
possibility by numerical diagonalization. 
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